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We analyze the conditions that enable acceleration of particles to ultra-high energies, ~ 10^" eV 
(UHECRs). We show that broad band photon data recently provided by WMAP, ISOCAM, Swift 
and Fermi satellites, yield constraints on the ability of active galactic nuclei (AGN) to produce 
UHECRs. The high energy (MeV - GeV) photons are produced by Compton scattering of the 
emitted low energy photons and the cosmic microwave background or extra-galactic background 
light. The ratio of the luminosities at high and low photon energies can therefore be used as a probe 
of the physical conditions in the acceleration site. We find that existing data excludes core regions 
of nearby radio-loud AGN as possible acceleration sites of UHECR protons. However, we show that 
giant radio lobes are not excluded. We apply our method to Cen A, and show that acceleration of 
protons to ~ lO'^" eV can only occur at distances > 100 kpc from the core. 



PACS numbers: 98.54.Cm, 95.85.Ry, 98.70.Sa 

The origin of ultra- high energy cosmic rays (UHECRs) , 
with energies > 10^^'^ eV, is still under debate. Astro- 
physical sources of UHECRs are limited by two require- 
ments: a strong magnetic field is needed to confine the 
accelerated cosmic rays, while the magnetic field cannot 
be too strong in order to avoid excessive synchrotron ra- 
diation and photo- meson energy losses In addition, 
energy losses by photo-meson production through scat- 
tering off the cosmic microwave background (CMB) limit 
the distance of UHECR sources to < 100 Mpc (the so- 
called GZK cutoff) 0, i]. Several possible sources of 
UHECRs that fulfill these constraints are discussed in 
the literature. The leading candidates are gamma-ray 
bursts (GRBs) [#|6|, low luminosity gamma-ray bursts 
and hypernovae [THlOj , and active galactic nuclei ( AGNs) 
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The possibility that AGNs are the main source of UHE- 
CRs is challenged by the need for an unusually high pho- 
ton luminosity. The required strong magnetic field, com- 
bined with the assumption that the energy density in the 
photon field is at least comparable to the energy density 
in the magnetic field, constrain the power output from 
the source to be L > lO^^F^/J-i ergs"! Here, 
F and /3c are the Lorentz factor and characteristic veloc- 
ity within the source. Only very few AGNs within the 
GZK horizon fulfill this requirement. A detailed analysis 
of those AGNs whose position was found to be correlated 
with the arrival direction of UHECRs ^] , resulted in a 
similar conclusion [21j. 

The underlying assumption of the above mentioned 
constraint is that the energy density in the photon field 
is comparable or higher than the energy density in the 
magnetic field. It is indeed difficult to obtain direct con- 
straints on the value of the magnetic fields in AGNs and 
test this assumption. 

In this paper, we re-analyze the constraints on AGNs as 
UHECR sources by combining low energy (radio-optical) 



and high energy MeV band) data, enabled by recent 
Swift and Fermi observatories. These data provide con- 
straints on the magnetic field strength, and avoid the 
need to adopt the equipartition assumption. In addition, 
the broad band data itself places constraints on the al- 
lowed region within the AGN that can produce UHECRs. 
The constraints are based on observable quantities, and 
thus overcome the inherent uncertainty in estimating the 
magnetic field strength. 

Basic requirements. Assuming that UHECRs accel- 
eration results from electromagnetic processes within an 
expanding plasma, the requirement that the acceler- 
ated particles are confined to the acceleration region is 
equivalent to the requirement that the acceleration time, 
tacc — TjE"'' /{rZeBc)\E^, is shorter than the dynamical 
time, tdyn — r/Tj3c. This condition implies 
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Here, is the observed energy of the particle, Ze is its 
charge, B is the magnetic field and r is the characteristic 
size of the acceleration region. 77 > 1 is a dimensionless 
factor, whose exact value is determined by the uncertain 
details of the acceleration mechanism |56j. Here and be- 
low, we use the convention Q = IQ^Qx in cgs units. 

A second condition is that the acceleration time 
is shorter than all relevant energy loss time scales. 
Energetic particles can lose their energy via syn- 
chrotron emission, on a time scale icooi.syn — 
{6TrmyTA^)/{aTmlB^E°''Z''). Here, ar is Thomson's 
cross section, rup and rrie are the proton and electron 
masses and Arrip is the mass of the nucleon (for iron 
nuclei, ^ = 56 and Z = 26). The requirement face < 
^cool, syn results in 
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Thus, both protons and iron nuclei can be accelerated to 
the highest observed energies, ^ 10^° eV, provided that 
the strength of the magnetic field at the acceleration site 
does not exceed B < few - few tens G. 

In addition to synchrotron energy losses, energetic par- 
ticles can in principle lose their energy by interacting 
with the ambient photon field and with other nuclei. 
Interaction with the photon field can result in energy 
losses by Compton scattering (which is typically a negligi- 
ble energy loss channel for UHECRs) , photopair produc- 
tion (Bethe-Heitler process), photo-production of mesons 
(mainly pions), and photodisintegration of heavy nuclei. 
Interaction with other nuclei can lead to spallation of 
heavy nuclei. These processes were studied in details 
for radio-loud AGNs in |22| , and for radio-quiet AGNs in 
[H]. The resuhs of both these works show that - lO^" eV 
particles will survive all energy losses, provided that the 
acceleration site is located at ~ few - few tens of parsecs 
from the core. 

Assuming that electrons are being accelerated at the 
same acceleration site of the UHECRs, synchrotron emis- 
sion is expected. Indeed, in many AGNs the peak of the 
synchrotron emission is clearly identified at the radio or 
infrared bands. The peak of the synchrotron energy flux 
can be approximated by assuming that the energetic elec- 
trons have a characteristic Lorentz factor 7e, 



one obtains 
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Here, ne is the number density of energetic electrons, V 
is the volume of the emitting region, is the luminosity 
distance, ctt is the Thomson cross section and V = [r(l — 
cos(d°'')]~^ is the Doppler factor for an observing angle 
9°^ [57|. For an observer within the light cone, 6°'' < 
max(r~^, 0,et)i P — r, where 9jet is the physical jet 
opening angle |58|. 

Recent Swift and Fermi observations revealed a sec- 
ond peak in the S2ectra of many AGNs, centered at 
~ MeV [e.g., |2J-|26[. This second peak results from 
inverse- Compton (IC) scattering of low energy photons 
by the energetic electrons; a hadronic origin can be ex- 
cluded by the lack of a neutrino counterpart [27[. Let 
us denote the flux ratio of the IC and the incoming 
photons peak energies by ^ = (j'Fi.)poak,ic/(i^-Pi^)poak,in. 
We do not assume that the incoming photons are nec- 
essarily the synchrotron photons, and we thus denote 
their frequency by fin, which can generally be different 
from Vsyn- In an IC process, the outgoing photon en- 
ergy is Vie = (4/3)7g^'m, and the outgoing monochro- 
matic (number) flux ratio of IC to synchrotron emission 
is (Fi.)peak,ic/(-F'i/)pcak,in = T ~ MueCTT ■ Approximat- 
ing the volume of the acceleration region as cylindrical, 
V = Trr^A/, and substituting these results in equation 
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Combined with the constraint set in equation ([1}, one 
flnds 



(5) 



Emission from AGN cores. A primary source of emis- 
sion in AGNs is the core, or inner parts of the jet. 
Due to the high brightness temperature of the ra- 
dio emission from this region, the seed photons to IC 
scattering are the synchrotron photons. Thus, R = 

(l'-Fi/)poakJc/(j^-F^)pcak,syn = Y ■, whcrC Y = (4/3)7eT is 

Compton parameter for an optical depth t. Equation 
([5|) can be written as a constraint on the minimum syn- 
chrotron luminosity that is needed for a source to be able 
to produce UHECRs, 
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Figure [1] shows the allowed region of the synchrotron lu- 
minosity as a function of Y (59|. Two nearby sources, 
Cen A and M87, have good enough spectral coverage for 
a reliable measurement of the synchrotron peak and the 
Y parameter. Both fall within the excluded region for 
acceleration of UHECRs to 10^° eV, unless UHECRs are 
composed of heavy nuclei. We point that for F ^ 1, 
the electrons can rapidly cool, leading to a broad band 
emission spectrum [28j. In this scenario, it would not 
be possible to identify both the synchrotron and the IC 
peaks, and the analysis could not be carried out. 

Emission from AGN lobes. While acceleration of 
UHECRs in the inner regions of AGNs are excluded by 
the data on Cen A and M87, jets from radio loud AGNs 
extend out to hundreds of kpc. Thus, additional possible 
acceleration sites exist in the turbulent outflow or at the 
termination shock of the giant AGN lobes. 

Due to the large distance from the source, the seed 
photons for the IC scattering in this region are not 
the synchrotron photons, but rather photons originat- 
ing from the cosmic microwave background (CMB) or 
extra- galactic background light (EBL). The incoming 
flux of the seed photons is therefore (j^Fiy)pcak,in = 
7rr^cUex/47r(i^, where r = i?acc is the characteristic size 
of the acceleration region, and Uox is the energy density 
of the CMB/EBL radiation fields. Substituting this re- 
sult in the definition of R in equation ([5]), one obtains the 
minimum radius which allows acceleration of UHECRs, 
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FIG. 1: Minimum peak luminosity of synchrotron emission in 
a source that is capable of accelerating UHECRs to 10^" eV, 
as a function oi Y = (!^-Fi,)peak,ic/(j^^l)peak,syn- The shaded 
(light yellow) region to the left of the solid line is excluded for 
{r)E!il/Z)pV = 1. The dashed line is for iron nuclei (Z=26). 



The measured values for the FR-I radio 
M87 are marked, based on data in Refs. 



galaxies Cen A and 

MM- 



In evaluating the minimuni radius in the second hne 
of equation (O, we conservatively used ucmb — 4 x 
10~^^ ergcm"'^. While several models exist for the en- 
ergy density of the EBL light [2^, it is generally smaller, 
UEBL ^ 10"^'' ergcm"^. Thus, if EBL photons are the 
target photons, a tighter constraint on the acceleration 
radius is obtained. The allowed region is plotted in Fig- 
ure [1160] 

Acceleration of UHECRs in Cen A. At distance of 
« 3.7 Mpc [U, Cen A is the nearest AGN. Due to 



its proximity, it was long been considered as a possible 
source of UHECRs [32, 33]. A renewed interest in its 
ability to produce UHECRs was prompted recently by 
the realization that the arrival directions of two and po- 
tentially even four UHECRs coincide within errors with 
its position 13- 3|. By analyzing VLA data, (s^l con- 
cluded that the inner jet is trans-relativistic, /3 ~ 0.5, 
with Doppler shift V > 1. A similar conclusion was re- 
cently drawn by the Fermi team [25[ , where V few was 
estimated, based on broad-band fits to the spectrum. 

Emission from Cen A's giant radio lobes, which ex- 
tend to > 300 kpc, was recently resolved by the Wilkin- 
son Microwave Anisotropy Probe (WMAP) [s^. By 
analyzing the data and estimating the magnetic field 
strength and the jet magnetic luminosity, it was sug- 
gested by Ref. 138] that the giant lobes are possible 
acceleration sites of UHECRs. The giant lobes were 
later resolved by the Fermi observatory [30| , constraining 

(i^-Fi/)peak,Ic/(i^-P'i/)poak,syn > 10. 

The results presented in Figures [T] and [51 which are 



FIG. 2: Minimum radius of acceleration region that allow 
acceleration of UHECRs to 10^" eV, as a function of the 
ratio (!^Fi,)pcak,ic/('^-F'i')peak,syn. The shaded (light yellow) 
region below the solid line is excluded for {riE2o/Z)f3V — 
1. The dashed line is for iron nuclei (Z=26). The 
dash-dotted (vertical) line gives the minimum ratio of 
(!^-F,.)pcak,ic/(j^fl)peak,syn > 10, measured for Cen A Q. 
Thus, Cen A can be the source of UHECRs, provided that 
the acceleration occurs beyond ~ 110 kpc. 



based on broad band data, show that the core region 
of Cen A does not fulfill the requirements that would 
enable it to accelerate protons to 10^" eV. However, the 
results presented in equation ^ and Figure [2] do not rule 
out the proposal [sl] that acceleration to 10^'^ eV could 
occur in the outer regions of the giant lobes. Using the 
additional data enabled by Fermi, we add a constraint 
on the minimum acceleration radius to be > 110 kpc. 

An upper limit on the flux of UHECRs from Cen A is 
obtained by assuming that all 13 events seen within 18° 
of Cen A [39| indeed originate from its lobes. This implies 
an observed cosmic rays flux J' « 2 x 10^^^ m^^s^^sr^^, 
which is translated into luminosity ~ ATid\jE°^ — 
-[^q38.5 ei-gg-i^ This value is about three orders of magni- 
tude less than the synchrotron luminosity of Cen A (see 
Figure [T|) , implying energetic consistency, even when in- 
terpolating the UHECRs flux to lower energies, provided 
that the power law index is not significantly greater than 
2. 

Discussion. We derived constraints on the ability of 
AGNs to accelerate UHECRs. These constraints are set 
as limits on observable quantities, in particular the syn- 
chrotron peak luminosity and the ratio of IC to syn- 
chrotron peaks. We showed in equation ([5]) and Figure 
[T]that the inner core regions of Cen A and M87 do not 
fulfill the necessary requirement, and therefore cannot 
be sources of UHECRs at 10^° eV, unless UHECRs are 
composed of heavy nuclei. 

On the other hand, we showed that outer regions of gi- 



4 



ant radio lobes, which are frequently seen in radio galax- 
ies, are possible acceleration sites. Using WMAP and 
Fermi data, we derived in equation ([7]) and Figure [5] the 
constraint that the acceleration site of the closest radio 
galaxy, Cen A, must be > 110 kpc. 

A common criticism of the possibility that AGNs might 
be the main sources of UHECRs is the lack of sources 
with bolometric luminosity Lboi > lO'^^ergs"^ within 
the GZK horizon [e.g., [io. Ell and references therein], 
requiring potentially bright flares [i^ [but see li^j • How- 
ever, the origin of this requirement lies in the assumption 
of equipartition of energy between photons and magnetic 
field, and is therefore not of general validity. 

Indeed, due to the difficulty in measuring the mag- 
netic field stren gth equipartition (or close to it) is often 
. Based on this assumption it was 
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assumed 1131. 1 161 . 
claimed by [16| that AGN cores are plausible acceleration 
sites. However, this claim did not consider the limitation 
of the observed luminosity in nearby AGNs. Moreover, as 
was shown in 13[ , energetic protons cannot escape inner 
regions without significantly energy loss due to photo- 
pion production, which limits the maximum energy of 
observed UHECRs. It was further pointed out by Ref. 
(iil that acceleration of CRs to high energies during flar- 
ing activities is limited due to insufficient residence time 
in the accelerated region. 

We presented a refined analysis, based on identifica- 
tion of the synchrotron and the IC peaks in broad band 
spectra of AGNs, which became available for Cen A and 
M87 through observations by the Swift and Fermi ob- 
servatories. By doing so, we removed any dependence on 
the limited validity of the equipartition assumption. This 
analysis applies to 'dark' sources, namely those sources 
in which the energy density in the magnetic field is much 
larger than the energy density in the photon field. This is 
indeed the case for the giant radio lobes of Cen A, where 
at ^ 100 kpc, the data shows that the ratio of magnetic 
to photon energy densities is ^ 10^ [38]. 

A detailed spectral analysis of many blazars was car- 
ried by 47|. In this analysis, a leptonic model (syn- 



chrotron and synchrotron self-Compton) was used to fit 
the broad band spectra of 73 nearby blazars (all of which 
are outside the GZK horizon). The obtained fits thus do 
not rely on the equipartition assumption. This analysis 
is therefore very similar to the one presented here, even 
more detailed. The fact that many blazars were found to 
have similar spectral shapes thus hints towards a possible 
generalization of our key results to many nearby AGNs. 

Using values obtained in a similar fitting to the broad- 
band spectra in Cen A 48|, it was claimed by [i^ that the 
core region of Cen A is excluded as an acceleration site 
of UHECRs. This conclusion is similar to the conclusion 
derived here, based on a much newer and better quality 
Fermi data. The analysis carried by Ref. [i^, however, 
did not consider the external giant radio lobes. 

In the analysis presented here we use a simplified as- 



sumption, namely that the electrons are mono-energetic. 
Thus, we ignore the energy distribution, which is needed 
for full spectral fits. Nonetheless, including a broad band 
spectral fitting, as was done, e.g., in Ref. [Hj does not 
change any of our conclusions. This is because only 
a small fraction of the electrons are being accelerated 
to high energies, and most of the emission occurs near 
the electron energy responsible for the synchrotron peak. 
This was verified by comparing numerical modelingof 
GRB spectra to simplified analytical approximations [50| . 

While our analysis excluded the core regions of Cen A 
and M87 (and plausibly other nearby AGNs) from being 
sources of UHECRs, giant fiares may change this conclu- 
sion. However, no measurement of the Compton param- 
eter during flares exists. If such a measurement becomes 
available, it could further constrain the ability of AGN 
cores to produce UHECRs during flaring activities. 

While existing data excludes the cores of Cen A and 
M87 from being potential acceleration sites of protons to 
10^'' cV, this is not the case if UHECRs are composed 
of heavy nuclei. Recent data collected by the Pierre 
Auger Observatory (PAO) show clues for heavy composi- 
tion [sH. This data, however, is in contrast to the results 
obtained by the HiRes detector js^l , and is therefore ten- 
tative. We point out though, that if indeed UHECRs 
are composed of heavy nuclei, not only the cores of radio 
loud AGNs are potential acceleration sites, but so many 
other objects, such as radio quiet AGNs [23|. 
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